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The four-body problem for an interacting two-species Fermi gas is solved analytically in a confined 
quasi-one-dimensional geometry, where the two-body atom-atom scattering length a aa displays a 
confinement-induced resonance. We compute the dimer-dimer scattering length add, and show that 
this quantity completely determines the many-body solution of the associated BEC-BCS crossover 
phenomenon in terms of bosonic dimers. 
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Cold atomic quantum gases continue to attract a lot 
of attention due to their high degree of control, tunabil- 
ity, and versatility. A main topic of interest has been 
the exploration of the BEC-BCS crossover in fermionic 
systems 0, & & IE El ■ I n two or three dimensions, this 
is still a controversial and not completely settled issue on 
the theory side 0,BEliEjj despite the qualitative agree- 
ment between mean-field theories and experimental data. 
Notably, a similar (but different) crossover phenomenon 
has been predicted to occur in quasi-one-dimensional 
(ID) systems [Hill, where a cylindrical trap leads to 
a confinement-induced resonance (CIR) [H Hi in the 
atom-atom interaction strength, analogous to the mag- 
netically tuned Feshbach resonance [8|. In contrast to 
what happens in 3D, one always has a two-body bound 
state ('dimer') in ID, regardless of the sign of the 3D 
atom-atom scattering length a. We solve the fermionic 
four-body problem in the confined geometry, and com- 
pute the dimer-dimer scattering length add throughout 
the full BCS-BEC crossover, on each side of the CIR. 
On the 'BEC side, we establish contact to results for 
the unconfined case ^5|, while on the 'BCS' side, a sim- 
ple Bethe Ansatz calculation provides exact results. The 
three-body problem has no trimer solution [l6l | , and thus 
the full many-body crossover solution can be expressed 
in terms of add alone and is thereby solved completely in 
this Letter. Since ID atomic gases can be prepared and 
probed thanks to recent advances our predic- 

tions could be observed in state-of-the-art experiments. 

We assume two fermion hyperfine components (de- 
noted by til) with identical particle numbers N-\ = 
Ni = N/2, interacting only via s-wave interactions. 
At low energies, the pseudopotential approximation [2(j 
for the 3D interaction among unlike fermions applies, 
V(r) = (4Trh 2 a/m )5(r)d r (r-) (m is the mass). We 
consider the transverse confinement potential U c (r) = 
moU!^_(x 2 +y 2 )/2, with lengthscale a± = (2h/m ujj_) 1 ^ 2 . 
The solution of the two-body problem 0, Q| reveals 
that a single dimer (composite boson) state exists for ar- 
bitrary a, where the dimensionlcss binding energy Qb 



and (longitudinal) size as, 
E B 



2hw 



K/2a B ) 2 > 0, 



(1) 



follow from ^(1/2, Qb) = — aj_/a with the Hurvitz zeta 
function. For an experimental verification, see Ref. |l8j |. 
For a±/a — > — oo, the BCS limit with £1^ ~ (a/a±) 2 <C 1 
and as — Oj_/2|a| is reached, while for a±/a — » +oo, the 
dimer (or BEC) limit emerges, with fls — (a±/2a) 2 ^> 1 
and as — a. The atom-atom scattering length is 

a aa = a x (C-a x /a)/2, C = -((1/2) ~ 1.4603. (2) 

For low energies, this result is reproduced by the ID 
atom-atom interaction V aa (z,z') = g aa S{z — z') with 
g aa = -2h 2 /m a aa jH]. The CIR (where g aa ±oo) 
takes place for aj_/a = C, which is equivalent to Qb = 1- 
In this paper, we solve the ID fermionic four-body (Till) 
problem and show that this also solves the A^-body prob- 
lem for arbitrary Qb in the low-energy regime. 

Let us first discuss general symmetries of the four-body 
problem. We denote the positions of the t (I) fermions by 
xi.4 (x.2,3), respectively, and then form distance vectors 
between unlike fermions, ri = xi — X2, r2 = X4 — X3, and 
r + = xi — X3, r = X4 — X2. The distance vector between 
dimers {12} and {34} is R./V2 = (x x + x 2 - x 3 - x 4 )/2. 
After an orthogonal transformation, the center-of-mass 
coordinate decouples and the four-body wavefunction \f 
depends only on ri 2 and R. With respect to dimer in- 
terchange, ^ is symmetric, 



*(ri,r 2> R) = *(r 2 ,ri,-R) J 
while under the exchange of identical fermions, 

*(r l5 r 2 ,R) = -*(r±,r T ,±(n -r 2 )/V2). 
The four-body Schrodinger equation then reads 



^ I A 

A P1 

m 



A r2 +A R ) + (7 c (r 1 ) + C/ c (r 2 ) 
U C (R)+V(r 2 ) -e]* = - V(Ti)V. 



(3) 
(4) 

(5) 



i=i. ± 
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The pscudopotcntials on the r.h.s. are incorporated via 
Bethe-Peierls boundary conditions imposed when a dimer 
is contracted, e.g., 

^(n.ra.RJI^o^^^a-ri/o). (6) 
47rri 

All other boundary conditions can also be expressed in 
terms of /(r,R) using Eqs. © and (0J, where 



/(r,R)=/(-r,-R), 



(T) 



expresses (parity) invariance of Eq. (JSJ) under r^a — > 
— ri 2 and R — > R in combination with Eq. ©. In 
order to appreciate the importance of Eq. Q, it is in- 
structive to expand /(r, R) in terms of the single-particle 
eigenfunctions tp\(r) and the two-body scattering states 
<&a(i") in the presence of the confinement, 



/(r,R) = ^/^$ /1 (r)^(R). 



(8) 



The quantum numbers A include the ID momentum k 
|2lj . the (integer) angular momentum m, and the radial 
quantum number n = 0,1,2, .. .. Explicit expressions for 
if)\ and <&a can be found in Refs. |13l 13. While both 
have the same energy E\, the $a now include the dimer 
bound state (denoted by A = 0) $o( r )- For relative longi- 
tudinal momentum k of the two dimcrs, the total energy 
is (excluding zero-point and center-of-mass motion) 



E 



h 2 k 2 

2mo 



(9) 



We consider the low-energy regime ka± < 1, where the 
relative dimer motion is in the lowest transverse state 
(n = m = 0) when dimers are far apart. We then have 
to deal with a ID dimer-dimer scattering problem in this 
'open' channel, where the asymptotic ID scattering state 
fo(Z) for \Z\ 3> max(ai, |a aa |) follows from Eq. ijSJ) as 

/(r, R) = $o(r)Vx,oo (VX 2 +Y 2 ^) fo(Z), (10) 

where ip± t oo is the transverse part of ipn=o.m=o- The 
symmetry relation (JTJ) now enforces fo(Z) = fo(—Z), re- 
flecting the fact that two (composite) bosons collide, i.e., 

fo(Z) = e-' lk W + (1 + 2/(/c))e ife|z| . (11) 

As long as only s-wave scattering is important, symme- 
try considerations thus rule out odd-parity solutions nor- 
mally present in ID scattering problems [l3L Il6j. This 
crucial observation implies that, assuming analyticity, 
the ID scattering amplitude can be expanded in terms of 
a ID dimer-dimer scattering length add H3- 



/>) = -l + ika dd + 0(k 2 ). 



(12) 



For |/cadd| <C 1, this also follows from the zero-range ID 
dimer-dimer potential 



V dd (Z,Z')=g dd 5(Z-Z'), 



2h 2 



9dd 



(2m )a dd 



(13) 



We stress that Eq. (|13|l holds for arbitrary aj_/a, and 
therefore ID dimer-dimer scattering at low energies is 
always characterized by a simple <5-interaction. 

Let us then analyze the BCS limit, fls <C 1, where 
the scattering problem is kinematically ID on length- 
scales exceeding a±. Projecting Eq. © onto the trans- 
verse ground state, the ID Schrodinger equation for 
four attractively interacting fermions reads with a aa = 
a 2 ± /2\a\ > a±, see Eq. ©, 



2m E 



i=l 



— * = 0, (14) 



where the second sum excludes identical fermion pairs, 
corresponding to {14} and {23}. The Bcthc Ansatz 
expresses the wavefunction as a sum of products of plane 
waves |23j . Let us choose the momenta a aa ki^ = ^fi — 
u/2 and a aa k3^ = T* + u/2 to describe dimer-dimer scat- 
tering, and measure lengths in units of a aa . The energy 
of this state is E = h 2 (—2 + u 2 /2)/(moa 2 a ) and u the 
relative momentum of the two dimers. Up to an overall 
normalization constant, the wavefunction in the domain 
T>i = {(zx, Z4) < (23, Z2)} must then be given by 



_ „-(z2+Z3-Z4-zi) I iu(z 2 +z 4 -z 3 -z 1 )/2 



\J> X = g — (*2-t-Z3— Z4 — 2l) /g- 

_ e iu(z 2 +z 1 -z 3 -z 4 )/2 __|_ e iu(z 3 +z 1 -z 2 -z 4 )/2 



iu{z 3 +z 4 -z 2 -zx)/2 



to ensure a normalizable and antisymmetric solution un- 
der exchange of identical fermions. Consider next a sec- 
ond domain, 2? 2 = {z\ < Z3 < Z4 < z%}, where 2:3 and 
Z4 are exchanged compared to T>\ . At the boundary be- 
tween T>\ and T>2, ^3 = Z4, Eq. 114|) implies "fi = \&2 and 



(^ 3 -a X4 )(* 1 -* 2 ) = -4* 1 



* 2 = 2Rc 



-(z 2 +z 3 —z 4 -z 1 ) _ 



leading to 



iu(z 2 +z 4 -z 3 -z 1 )/2 



2 + ill 



/ 2 

_|_ e ~(z 2 +Z4 — 23 — 21) * e iu(z 2 +z 3 — 24 — 2i)/2 

\2 + iu 

_ 2g iu ( Z 2+2l-2 4 -23)/2) 



The wavefunction in other domains can be found in a sim- 
ilar manner. As a result, for a large dimer-dimer distance 
Z, ^ oc e _|z+l e" |2 - | /o(^) J where e" |z±l is the ID wave- 
function of the dimer {13} and {24}, respectively. This 
result shows explicitly that even in the BCS limit, the 
two dimers are not broken in the collision even for large 
k. There is no coupling to additional fermionic states, 
and the composite nature of the dimer is not apparent in 
<I>. The ID scattering state fo(Z), see Eq. (fTTfl - has the 
exact scattering amplitude 



1 



1 + ika dd 



add 



a ± 
4|o|' 



(15) 



which reproduces the full scattering amplitude derived 
from Eq. (|13f) and not just the first order as in Eq. i|12[) . 
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The bound state at imaginary k predicted by Eq. i|15|) 
is however unphysical, since the corresponding Bethe 
Ansatz solutions are then not normalizable. It would cor- 
respond to a non-existent bound four-fermion (tetramer) 
state, and hence Eq. (|15fl is restricted to the real axis. 

Let us now turn to the many-body problem, starting 
with the BCS limit. Since dimers are not broken in the 
collision, the ground state can be described in terms of 
N/2 bosons ('bosonization') with the interaction (|13|l and 
a dd = a aa/2- The attractively interacting Bose gas is sta- 
bilized by the real-Zc restriction, implying the omission of 
many-body bosonic bound states. Bosonization is pos- 
sible for pa± < 1, since typical momenta are k f=a p for 
total ID fcrmionic density p. This reasoning immediately 
leads to the famous Lieb-Liniger (LL) equations j2f| , 



E ° ft o ^ 

— = -ftUj_llB + ~ 

N p 

A'o 



K 



2nf(k) = 1 



dp ■ 



n 2 k 2 

m 



ldJ-K 4/a dd + (p-k) 2 '' 



(16a) 



(16b) 



where Eq is the ground state energy and Kq is fixed by 
p/2 = dkf(k). Notably, since add = a aa /2, the LL 
equations coincide with Yang-Gaudin equations for N at- 
tractively interacting ID fermions, thereby exp laining a 
deep connection noticed previously 0, 12, 261. Moving 
towards the dimer limit, Eq. (|13|) still applies, but now 
only for sufficiently small k such that Eq. (|12|) holds, and 
o-dd 7^ a aa/2. For add ^ a±, one leaves the BCS regime 
and enters the 'crossover regime', while (once add < 0) 
the dimer regime is realized for \add\ a±. Within the 
crossover regime, \add\ ^ cl±, we have hard-core bosons 
that can effectively be fermionized |lll ll2T | , again imply- 
ing typical momenta k p. For pa± < 1, the condition 
|fccw| ^ 1 imposed by Eq. I|12|) is therefore safely ful- 
filled throughout the crossover regime. Finally, in the 
dimer limit, a < aj_, fermions form very tightly bound 
dimers. The confinement can then not influence the four- 
body collision, which is therefore described by a 3D zero- 
range interaction with a dd w 0.6a |l5l |. However, for 
dimer-dimer distance larger than «i, dimers eventually 
must occupy the transverse ground state, see Eq. (|TU|> . 
In effect, for pa± < 1, we recover a ID (bosonic) two- 
body problem, where Eq. @ gives the answer (exact for 

n B > i), 



add 



"rcd,_L 

2(0.6a)' 



a Ied ,± = (h 2 /m Lu ± A 2 , (17) 



where a ro d.± is the transverse oscillator length for dimers. 
To summarize this discussion, we have shown that (a) as 
long as the single condition ka± < 1 holds, dimer-dimer 
scattering is described by Eq. I|13|) for arbitrary a±/a, 
and (b) knowledge of add and hence the solution of the 
ID four-body problem is sufficient to completely solve 
the ID BEC-BCS many-body problem for dilute systems, 
pa^ < 1, in terms of the LL equations l|16|) . 

Next we discuss the ID four-body problem. Enforcing 
the boundary condition © or the other equivalent ones, 




FIG. 1: Scattering length add as a function of SIb- Dashed 
curves give exact limiting results, the solid curve interpo- 
lates by adding these. Inset: Same but neglecting all closed- 
channel excitations. Here the solid curve gives the exact re- 
sult. 



Eq. (JSJ leads to an integral equation for /(r, R) [T3.ll^. 
Using Eq. ©, some algebra [23 yields 



-C(i/2,n B ) 



Uu (is) 



47r7t 2 a_L \- rll ' v ' , 

^'=E / drdRG B ^_ Eu ((r±V2R)/2,0) 
r 

3, 




e( T ^]^(r)^(R). 



The two-body Green's function Ge^, 0) can be found 
in Ref. ^(|. The two degrees of freedom in /(r, R) im- 
ply two different types of 'closed' channels that may be 
excited in a dimer-dimer collision: (i) scattering states 
above the bound state for each dimer [corresponding to 
r or p in Eq. JHJl], and (ii) excited states in the trans- 
verse direction for the relative motion of two dimers [cor- 
responding to R or v in Eq. JHJl]. Neglecting both types 
of closed- channel excitations, Eq. (|18f) can be solved nu- 
merically for arbitrary a±/a as in Ref. [Tt| . The result 
is shown in the inset of Fig. 2] In addition, this approxi- 
mation allows to extract add hi both limits analytically: 
in the dimer limit, we find add = —KQa\/(2a) + 2n\a, 
where = 1/4 and n\ ~ 0.319, while in the BCS limit, 
a-dd = ^oal/M with ijq ~ 0.402. The exact (numeri- 
cal) result for arbitrary a±/a agrees to within ±0.05 in 
o-dd/a-i. with a simple interpolation formula obtained by 
simply adding these two limiting results. For practical 
purposes, the interpolation is therefore virtually exact. 
Let us then turn to the effects of closed-channel excita- 
tions. In the BCS limit, excitations of type (ii) are irrel- 
evant [l6l |. but type-(i) excitations are important. Their 
inclusion results in the exact value rj = 1/4, see Eq. <|15ll . 
which also follows from the solution of Eq. i|18[l including 
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FIG. 2: Squared ratio of breathing and dipole mode fre- 
quency as a function of —ax/a. Here we have chosen 
Nuj z /lj ± = 1/3. 

typc-(i) excitations |27^ . In the dimer limit, inclusion of 
the closed channels leads to the correct value kq r* 0.83, 
see Eq. <|17[) . Incidentally, the two excitation types can 
be disentangled [23], and we find a d ® ~ 0.66a by just 
neglecting type-(i) excitations, which is already close to 
the exact value a dd « 0.6a [15|. Type-(ii) excitations 
are obviously important in the dimer limit, which may 
be valuable input for diagrammatics (TJ, |23| • The exact 
limiting results for add are shown in the main part of 
Fig. ^ as dashed curves. For the full crossover, the addi- 



tive interpolation formula is again expected to be highly 
accurate. Notably, this predicts add = for 0# sa 0.3. At 
this point, a CIR for dimer-dimer scattering occurs, see 
Eq. I|13|) , where the interaction strength gdd diverges and 
changes sign. Interestingly, the dimer-dimer CIR takes 
place at a different value for (and hence a±/a) than 
the atom-atom CIR. 

In experiments, quasi-lD regimes can be obtained in 
arrays of very elongated traps with a shallow confine- 
ment in the longitudinal direction. Typical trap frequen- 
cies are oj±/2it sa 70 kHz and lu z /2tt sa 250 Hz, with 
N sa 100 atoms per tube to ensure the ID condition 
N < uj ± /u; z pi]. The BCS-BEC crossover can be in- 
vestigated using a Feshbach resonance, which leads to 
changes in the density profile excitation gaps 
and ground state energy that can be probed via release 
energy p| and rf spectroscopy measurements 0] . A 
probably more precise approach is to measure collective 
axial modes. The dipole mode frequency is always u> z , 
irrespective of interactions. Using a sum rule approach 
[2^ . we calculated the frequency of the lowest compres- 
sional (breathing) mode from the mean-square size of the 

cloud bJ 2 = —2(d\n(z 2 )/dwl) , see Fig. [21 by solving 
Eqs. (|16(l using our results for add- Limiting values arc 
w = \/3u) z in the dimer limit, and ui = 2uj z both in the 
BCS limit and close to add = 0. Wc hope that this pre- 
diction will soon be tested. - This work was supported 
by the DFG-SFB TR12. 
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